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Abstract. Since the first experimental realization of Bose-Einstein 
condensation in cold atomic gases in 1995 there has been a surge of ac- 
tivity in this field. Ingenious experiments have allowed us to probe mat- 
ter close to zero temperature and reveal some of the fascinating effects 
quantum mechanics has bestowed on nature. It is a challenge for math- 
ematical physicists to understand these various phenomena from first 
principles, that is, starting from the underlying many-body Schrodinger 
equation. Recent progress in this direction concerns mainly equilibrium 
properties of dilute, cold quantum gases. We shall explain some of the 
results in this article, and describe the mathematics involved in under- 
standing these phenomena. Topics include the ground state energy and 
the free energy at positive temperature, the effect of interparticle in- 
teraction on the critical temperature for Bose-Einstein condensation, as 
well as the occurrence of superfluidity and quantized vortices in rapidly 
rotating gases. 



1. Introduction 

Bose-Einstein Condensation (BEC) was first experimentally realized in 
cold atomic gases in 1995 [2j [10]. In these experiments, a large number 
of (bosonic) atoms is confined to a trap and cooled to very low tempera- 
tures. Below a certain critical temperature condensation of a large fraction 
of particles into the same one-particle state occurs. 

These Bose-Einstein condensates display various interesting quantum phe- 
nomena, like superfluidity and the appearance of quantized vortices in rotat- 
ing traps, effective lower dimensional behavior in strongly elongated traps, 
etc. We refer to the review articles [9], El El [15] for an overview of the state- 
of-the-art of this subject and a list of references to the original literature. 

BEC was predicted by Einstein in 1924 [13] from considerations of the 
non- interacting Bose gas, extending the work of Bose [7] to massive parti- 
cles. The presence of particle interactions represents a major difficulty for a 
rigorous derivation of this phenomenon, however, as we shall discuss below. 

1.1. The Bose Gas: A Quantum Many-Body Problem. The quantum- 
mechanical description of the Bose gas is given in terms of its Hamiltonian. 
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For a gas of N bosons confined to a region A G R 3 , and interacting via a 
repulsive pair-interaction potential v, it is given by 



(1) 



The kinetic energy of the particles is described by the Laplacian A = V 2 , 
and we choose Dirichlet boundary conditions on dA for concreteness. Other 
boundary conditions could be used as well. The subscript i stands for the 
action in the ith particle coordinate x% G R 3 . Units are chosen such that 
h = 2m = 1, with m the mass of the particles. 

The Hamiltonian H acts on the Hilbert space of permutation-symmetric 
wave functions \I> G (£) N L 2 (R 3 ), as appropriate for bosons; i.e., square- 
integrable functions of N variables Xi G R 3 satisfying 

%1){XI,. . . ,X N ) =l/>(x n (!),...,X n ( N )) 

for any permutation tt of (1, 2, . . . , N). 

In the following, the interaction v will be assumed to be radial and non- 
negative. Moreover, it is sufficiently short range as to have a finite scattering 
length, which means that it is integrable outside some compact set. No other 
regularity assumptions will be made. In particular, v is allowed to have a 
hard core, which reduces the domain of definition of H to those functions 
that vanish whenever the distance between a pair of particle coordinates is 
smaller than the hard-sphere radius. A particular example of an interaction 
potential to keep in mind are pure hard spheres where, formally, v(x) = oo 
for | a; | < a, and v(x) = for \x\ > a. 

The present setup can be easily generalized to describe inhomogeneous 
systems in a trap. One simply adds trap potential 

N 

i=i 

to H, where V is a real- valued, locally bounded function with limiji^Qo V{x) = 
oo. The latter condition guarantees that the particles are confined to the 
trap, even in case A = R 3 . 

Similarly, rotating systems can be described adding the term 

N 



i=l 



to the Hamiltonian H, with Q G M 3 being the angular velocity and L = —ixA 
V the angular momentum operator. This term results from a transformation 
to the rotating frame of reference. 
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1.2. Quantities of Interest. In the following, we shall distinguish two 
types of questions that can be asked concerning of the behavior of Bose 
gases described by the Hamiltonian ([I]) above. 

• Thermodynamic quantities, like the ground state energy per unit 
volume, or the free energy density at positive temperature. Here one 
considers homogeneous systems and is interested in the thermody- 
namic limit N — ► oo, A — ► R 3 with the particle density g = N/\A\ 
fixed. 

Of particular interest is the notion of Bose-Einstein condensation, 
which concerns off-diagonal long-range order in the one-particle den- 
sity matrix (aJ (x)a(y)), and is expected to occur below a critical 
temperature. 

• Behavior of trapped systems in the ground state. One observes 
interesting quantum phenomena, like effective one-dimensional be- 
havior in strongly elongated traps, vortices in rotating systems, a 
bosonic analogue of the fractional quantum Hall effect in rapidly 
rotating gases, etc. 

Of particular relevance is the Gross-Pitaevskii scaling, where 
the ratio of the scattering length a to the diameter of the trap is 
OiN- 1 ). 

We shall discuss our current knowledge about answers to these questions, 
as far as mathematical physics is concerned, in the following sections. 



2. Homogeneous Systems in the Thermodynamic Limit 

2.1. The Ground State Energy of Homogeneous Bose Gases. Con- 
sider first the case of a homogeneous system in the absence of a trapping 
potential or rotation. The ground state energy density in the thermody- 
namic limit is given by 

e(o) = lim 7T7 fnf suecH (2) 

A— R3,JV/|A|— e l A l 

with H as in ([TJ. The existence of this thermodynamic limit is well under- 
stood for appropriate sequences of domains A approaching R 3 . See, e.g., 
Ruelle's book pg|. 

We will be particularly interested in the limit of low density, when the gas 
is dilute in the sense that a 3 g <C 1, where a denotes the scattering length of 
the interaction potential v. It is defined as 

4vra = inf I J (\ V(/>(|x|)| 2 + \v{x )(f>(\x\) 2 ^J dx : <p > , lhn <j>(r) = lj . 

(3) 

For bosons at low density, one expects that 

e(g) ss iirag 2 . (4) 
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This formula is suggested by considering the ground state energy of two 
bosons in a large region A, which is 87ra/|A|, as can be easily deduced from 
(|3|). Multiplying this by the number of pairs of bosons, N(N — l)/2, one 
arrives at (j4]). That this simple heuristics is correct is far from obvious, 
however. It fails for two-dimensional systems, for instance [35j[29]. 

The investigation of the ground state energy density e(g) goes back to 
Bogoliubov [6] in the 40s, and Lee, Huang and Yang in the 50s [19]. Dyson 
[TT] computed a rigorous upper bound that shows the correct leading order 
asymptotics (jlj) for hard spheres, but his lower bound was 14 times too small. 
His upper bound was later generalized to arbitrary repulsive interaction 
potentials in [24]. The correct lower bound was proved only in 1998 by Lieb 
and Yngvason [28]. We formulate this result as a theorem. 

Theorem 1 (Bosons at T = 0). As g -> 0, 



e(g) = Airag 2 + o(g 2 ) (5) 



Note that if one treats the interaction energy as a perturbation of the 
kinetic energy, naive perturbation theory would yield I J v instead of Aira. 
This is always too big, as ([3]) shows, and would even be infinite for hard 
spheres. In fact, the result ([5]) is non-perturbative in the sense that the 
scattering length a contains terms to arbitrary high order in the interaction 
potential v. 

It remains an open problem to establish the leading order correction to 
(jl]), which is expected to be given by the Lee-Huang- Yang formula [1~9] 



(128 \ 
1 + j^V^gj . (6) 

Recent progress in this direction was made in [TJJ and [27], where it was 
shown that © holds for certain density-dependent and appropriately scaled 
interaction potentials. The general question remains open, however. 

2.2. Homogeneous Bose Gas at Positive Temperature. At positive 
temperature T > 0, the appropriate quantity to consider is the free energy 
density, which is defined as 

f(g,T) = -T lim -LlnTrexp(-tf/T). (7) 

A— >M 3 , N/\A\—*g \ A\ 

For non-interacting bosons (i.e., v = 0), it can be calculated explicitly. We 
denote it by fo(g,T). It is given in terms of a Legendre transform as 



fo(Q,T) = sup 

At<0 



(8) 



Note that it has the scaling property 

/o(e,T) = ^Vo(i,T e - 2/3 ) 
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which follows from the fact that the only length scales in the problem are the 
mean particle spacing g^ 1 / 3 and the thermal wavelength T _1//2 , and hence 
fo depends, up to a prefactor, only on their ratio. 

From (JSj) it is easy to see that /o is not an analytic function of g (or T), 
and hence even a non-interacting Bose gas shows a phase transition. This 
transition is known as Bose-Einstein condensation, and occurs at a critical 
density 

/ <T> \ 3/2 

Sc(T) = C(l) 7Z) , 



where C, denotes the Riemann zeta-function. In fact, dfo(g,T)/dg = for 
Q > 6c- For g > g c (T), g — g c (T) is interpreted as the density of the Bose- 
Einstein condensate. 

For interacting gases, there are now three length scales to consider: the 
interaction range a, the mean particle distance £> _1//3 , and the thermal wave- 
length T^ 1 / 2 . For dilute systems, one considers the case 

a « 0-V3 ~ T-V 2 . 

In this regime, the free energy turns out to be the given by the following 
expression. 

Theorem 2 (Bosons at T > 0). For a 3 g C 1 we have 



f(g,T)=fo(g,T)+4ira(2g 2 -[g-g c (T)]l) + o(g 



0) 



where = max{t, 0} denotes the positive part. 



The lower bound in §§§ was proved in [3B]. An upper bound for smooth 
interacting potentials of rapid decay was later obtained in [3D], the more 
general case being still open. 

The error term in ([9J is uniform in T/g 2 ^ 3 for bounded T/g 2 / 3 , corre- 
sponding to the quantum regime. For T j g 2 ! 3 — > oo one obtains a classical 
gas, whereas for T/g 2 / 3 — > the system approaches the ground state. 

Note that for g < g c (T), the leading order correction compared to the 
ideal Bose gas is 8nag 2 instead of the ^irag 2 at zero temperature. The 
additional factor 2 is a result of the symmetry requirements of the wave 
functions and can be interpreted as an exchange term; this symmetrization 
applies only to particles outside the condensate, however, and this explains 
the subtraction of the square of the condensate density in ([9]). We also 
remark that without restricting to symmetric functions, the leading order 
correction compared with an ideal gas would be Airag 2 at any T > 0, just 
like at T = 0. 

The proof of Theorem [2] is long and technical and hence can not be repro- 
duced here. One of the key issues to understand is a certain separation of 
energy scales in the two terms on the right side of ([9]). In momentum space, 
these are 
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• large momenta \p\ ~ 1/a responsible for scattering of two particles 
at a distance ~ a of each other. 

• low momenta \p\ ~ T 1 / 2 <C 1/a, responsible for the thermal distri- 
bution distribution of the particles' kinetic energy. 

• Bose-Einstein condensation at momentum p = 0. 



2.3. Critical Temperature for BEC. As discussed above, the ideal, non- 
interacting Bose gas displays a phase transition above a critical density. 
Equivalently, BEC in the ideal gas occurs below the critical temperature 

T (n) - — ^ o 2 / 3 

c[6) ~ C(3/2) 2 /3^ • 

A useful characterization of BEC, applicable also for interacting systems, 
is in terms of the one-particle density matrix of the system. This density 
matrix is defined as 

T » 1 £± Slf Tfi'-%-*IT (10) 

where Tr^ -1 stands for the partial trace over A — 1 particle coordinates. 
Hence 7 is an operator on the one-particle space L 2 (M 3 ). Obviously 7 > 
and Tr7 = A, by definition. BEC is characterized by the fact that, in the 
thermodynamic limit, the integral kernel j(x, y) of 7 does not vanish as 
\x — y\ — * 00. This is also referred to as off-diagonal long range order. For 
non-interacting bosons, one can show that 

7 (f , y) = [g- g c (T)] + + £ 6 e-ri^VW 

n>0 ( 47m / T ) 3/2 

in the thermodynamic limit, with [t]+ = max{t, 0} denoting the positive 
part, and p, < the fi where the maximum in (JSj) is achieved. Hence the 
kernel j(x, y) has the following characteristics: 

• For T < T c (g), j(x,y) does not decay. In fact, limi£_.j7|_ +lx , 7(2;, y) = 
g — Q C (T), the condensate density 

• For T > T c (g), j(x,y) decays exponentially, like e~^~~^ x ~^ 

• For T = T c (g), 7(3;, y) decays algebraically. In fact, 7(2?, y) ~ \x — 

in this case. 

These features are expected to hold also for interacting Bose gases, al- 
though with a different value of the critical temperature T c (g). It is still an 
open problem to prove the existence of BEC for interacting gases, however. 
The only known case where BEC has been proved is the hard-core lattice 
gas at half filling [12], which is equivalent to the XY spin model. [3~T] 

Although there is no proof that T c ^ in the interacting case, an upper 
bound can be derived rigorously [39] : 
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Theorem 3 (Upper bound on T c ). For small a?g and some c > 0, 




where T, 



(0) 



C(3/2)2/3 



n 2/3 



i/ie critical temperature for the ideal Bose gas. 



More precisely, it is shown in [39] that 7(2;, y) decays exponentially if 

T > T c (0) (l + Cx /o^ 1/3 )- The proof uses a well-known Feynman-Kac repre- 
sentation of the partition function in terms of integrals over paths and sums 
over cycles in permutations. [16] 

There seems to be still no consensus in the physics literature concerning 
the correct power of the exponent of ag 1 ^ in the shift in critical temperature, 
or even the sign of c! Recent numerical simulations suggest that the shift 
should be linear in ag 1 ^ , with a positive c. This expected behavior of T c (g), 
as well as the upper bound of Theorem El are sketched in Figure CD 



T/p 2/3 



47T 

C(3/2) 2 /3 




Figure 1. The red line shows the rigorous upper bound 
on the critical temperature for BEC. The dashed line corre- 
sponds to the expected behavior based on numerical simula- 
tions. 



3. Trapped Bose Gases 

In the previous chapter we considered homogeneous Bose gases in the 
thermodynamic limit. Recent experiments with cold atoms consider inho- 
mogeneous gases in traps, however. That is, one can take A to be the whole 
of M 3 , but adds a trap potential Y^iLi to the Hamiltonian (pQ). A 

typical example, which describes the experimental situation rather well, is 
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a harmonic oscillator potential V(x) = lo 2 \x\ 2 , with cu > the trap fre- 
quency. More generally, the trapping frequencies in the three directions can 
be different, of course. 

A characteristic feature of these trapped gases is their response to ro- 
tation. One observes the appearance of quantized vortices [30} I14j. whose 
number increases with the rotation speed 

Even a rotating Bose gas can be described in a time-independent way, by 
going to the rotating reference frame. The only effect on the Hamiltonian 
is to add the term J2iLi ^ ' as discussed in the Introduction. To ensure 
stability of the system, the trap potential V has to increase fast enough at 
infinity to compensate for the centrifugal force in the rotating system. More 
precisely, we have to assume that 

Jim (V(x) - A x\ 2 ) = +oo . 

\x\— >oo V / 



3.1. The Gross-Pitaevskii Equation. The previous considerations sug- 
gest that dilute Bose gases close to zero temperature should be well described 
by the Gross-Pitaevskii (GP) energy functional [HI [33] 

£ GP [<f)] = U -A + V{x) -n-L 4)+ iirNa [ \cj)(x)\ 4 dx . (11) 
\ / Jm.3 

Its ground state energy is 



E GP (Na,n)= inf £ GP [(j)} 
\\4>h=i 



and any minimizer satisfies the GP equation 



A(p(x) + V{x)<j>(x) - Ct -L<f>(x) + 87rA^a|(/>(f)| 2 (/>(f) = n<t>{x) 



For a minimizer, N\<ft(x)\ 2 is interpreted as the particle density of the system. 
Hence the last term in (|lip is the natural generalization of the expression 
4tt(iq 2 to inhomogeneous systems. 

For $7^0 and axially symmetric V(x), the rotation symmetry can be 
broken due to the appearance of quantized vortices. More precisely, it was 
shown in [361 [37] that for all trap potentials V(x) that grow faster than 
quadratically at infinity, there exists a such that for all Na > the GP 
minimizers necessarily are not axially symmetric. In particular, there are 
many (in fact, uncountably many) GP minimizers! The symmetry breaking 
is due to the appearance of quantized vortices which can not be arranged in 
a symmetric way. Many interesting results have been obtained concerning 
the nature and distribution of these vortices in GP minimizers. We refer to 
P] and references therein. 
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3.2. Ground State Energy of Dilute Trapped Gases. In typical ex- 
periments on cold atomic gases, JV > 1, a C 1 (the length scale of the 
trapping potential V), but Na = 0{1). To get to this dilute regime, one 
writes 

v(x) = —^wix/a) (12) 
a z 

with w having scattering length 1. It is easy to see that v(x) then has 
scattering length a. The scattering length thus becomes a parameter, and 
we can write 

inf spec H = Eo(N,a,Q) . 

We note that the scaling (|12p is of course equivalent to a rescaling of the trap 
potential V while keeping v fixed. This latter procedure may seem physically 
more natural (as the trap potential is easier to adjust experimentally than 
the interaction potential) but we find it more convenient to fix V instead 
and scale v as in (|12j) instead. Our procedure corresponds to measuring all 
lengths in the system in units of the length scale of the trap potential. 

For dilute systems, one expects that Eo(N,a,Q) ~ NE GP (Na,&). The 
proof of this fact was given in [23] . 

Theorem 4 (Ground State Energy of Trapped Gases). For fixed g > 
and U G R 3 , 

wpT ^I (13) 

N^oo iv 

This theorem was previously proved in [24] for the case SI = 0. The main 
difficulty in the generalization to rotating systems comes from the fact that 
the permutation symmetry of the wave functions now becomes essential. 
While for non-rotating systems it is well known that the ground state for 
bosons coincides with the ground state without symmetry restrictions (as 
the latter is unique and positive, hence must be symmetric), this fact fails 
to hold for rotating systems. In fact one can show that (|13p fails to hold, in 
general, if the left side is replaced by the absolute ground state energy of H 
(viewed as an operator on L 2 (M 3Ar ), without symmetry restrictions). [37] 

3.3. BEC for Rotating Trapped Gases. In the previous subsection it 
was argued that the ground state energy of the GP functional (jlip is a good 
approximation to the ground state energy of H for dilute gases. For the 
corresponding ground state ^o(^i ; • • • >^a0> one would also expect that its 
one-particle reduced density matrix satisfies 

7o = Tr( JV - 1 )|* )(*o| ~ |0)<0| (14) 

with (j) a minimizer of the GP functional. (For convenience, the normaliza- 
tion of 70 has been chosen differently here than we did previously in (jlOp .) 
While this is indeed true in the non-rotating case [22], the rotating case is 
more complicated because of the non-uniqueness of the GP minimizers <ft. 
The best one can hope for is to replace the right side of (|14h by a convex 
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combination of rank-one projections onto GP minimizers. This is indeed the 
content of Theorem [5] below, which was proved in [23] . 

To state the following results precisely, it is necessary to introduce the con- 
cept of an approximate ground state. We will call a sequence of iV-particle 
density matrices (positive trace class operators on the iV-particle space with 
trace equal to one) an approximate ground state if their energy equals the 
ground state energy to leading order in N. Then we define the set T as 
the set of limit points of one-particle density matrices of such approximate 
ground states. More precisely, 

r = {7 : 3 sequence 7^, lim ^-Tr Hj N = E GP (g,n), Jam 7$ = 7 

N— >oo, Na— *g iv N^oo 

(15) 

where jff = Tr^ _1 ^7jv denotes the one-particle density matrix of 7at. 

Theorem 5 (BEC for Dilute Trapped Gases). The set T in ITR) has 

the following properties. 

(i) r C J\ is compact and convex. 

(ii) The extreme points r ext C T are given by GP minimizers, i.e., r ext = 
(10X01 : £ GP [<j>}=E GP (g,ti)}. 

(hi) For every 7 G T there exists a positive (regular Borel) measure d[x^, 
supported in r ex t with J r t <i/u 7 ((/>) = 1, such that 

7 = / rfM-y(^) 10X01 - (16) 

J Text 

Eq. (I16h is the natural generalization of (I14p to the case of multiple GP 
minimizers. It says that the one-particle density matrix of any approximate 
ground state is close (in trace class norm) to the convex combination of 
projections onto GP minimizers. 

Theorem [5] represents also a proof of the spontaneous breaking of the 
rotation symmetry in rotating Bose gases. An infinitesimal perturbation, 
e.g. of the trap potential V, leads to a unique GP minimizer and hence 
to 100% condensation, since the set T consists of only one element in this 
case. The quantized vortices are visible in the GP minimizer; they are a 
typical feature of superfluids. Theorem [5] can therefore also be interpreted 
as a proof of the superfluid behavior of rotating Bose gases. [25] 

3.4. Rapid Rotation. Consider now the special case of a harmonic trap- 
ping potential 

V{x) = \\x\ 2 . 

As discussed above, H is bounded below only for < 1. The results in 
the previous subsections are valid for fixed < 1. The question we would 
like to address in this final section is what happens as |0| — > 1? Denoting 



HOT TOPICS IN COLD GASES 



11 



eh = 0/| O | the unit vector in the direction of O, we can write 
- A + \\x\ 2 -Q-L= (-iV - A eh A x) 2 + ±|e n ■ x\ 2 +(e n - O) • X . (17) 



h 



The operator h has eigenvalues | , | , | , . . . , each of which is infinitely degen- 
erate. 

For low energies it makes sense to restrict the allowed wave functions to 
the kernel of h — | . This kernel is given by the Bargmann space [3] 

{f{z)e-W 2/4 , f : C -» C analytic} C L 2 {M 3 ) (18) 

where we identify the complex variable z with the plane perpendicular to f2. 
In particular, \x\ 2 = \z\ 2 + |eh ■ x\ 2 . Since the Gaussian factor is fixed, it is 
convenient to absorb it into the measure and think of the Hilbert space as 
a space of analytic functions only. It can easily be checked that the angular 
momentum operator eh ■ L acts on / as zd z . In particular, its eigenfunctions 
are z n , with eigenvalues n £ {0, 1,2,...}. 

We note that if one interprets eh as a homogeneous magnetic field, the 
Bargmann space (|18j) corresponds to the lowest Landau level in the perpen- 
dicular direction, multiplied by a fixed Gaussian in the longitudinal direc- 
tion. 

Now that we have identified the Bargmann space (|18p as the appropriate 
one-particle Hilbert space for rapidly rotating bosons, we have to come up 
with an effective Hamiltonian describing this system. The only term left 
in the one-particle energy (|17h . expect for a trivial factor |, is the angular 
momentum term (eh — ft) ■ L. If the range of the interaction potential is 
much shorter than the "magnetic length" 1, it makes sense to approximate 
the interaction potential by a ^-function, which becomes a bounded operator 
when projected to the Bargmann space. Writing the prefactor of the 6- 
function as 87ra, in accordance with previous considerations, we arrive at 
the effective Hamiltonian 



N 

H LLL := (1 - | d|) Y, + 8 ™ £ S ij ■ 

i=l l<i<j<N 



(19) 



It acts on the space of permutation-symmetric analytic functions f(z\, . . . , z^) 

which are square-integrable with respect to the measure 

We denote this space by B® N . The operator <$ 12 : B® 2 — ► B® 2 acts as 

(<W) (ZX,Z 2 ) = ^3/2 / + Z 2), \{Z\ + Zi)) 

which takes analytic functions into analytic functions. It is obtained by 
projecting 5{x\ — x%) onto B® 2 . 

Concerning the effective Hamiltonian (|19p . the following questions arise 
naturally: 
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(1) Can one derive H LLL rigorously from the full Hamiltonian H as 

— ► 1 and a — > 0? Such a rigorous derivation was indeed achieved 
in [21], where it was shown that if — > 1 for fixed N and fixed 
a/(l — |0|) (i.e., also a — > 0), then the ratio of the ground state 
energy of -ff, minus the trivial term |iV, to the ground state energy 
of iJ LLL goes to 1. Similarly, one can show that also eigenfunctions 
converge in the same limit. Uniformity in the particle number N is 
still an open problem, however. 

(2) What are the properties of H LLL , in particular concerning its spec- 
trum and corresponding eigenfunctions? Certain features of if LLL 
are expected to show some similarities to the fractional quantum 
Hall effect which occurs in fermionic systems. 

Concerning the latter question, let us first note that the two terms Cn 
and A 7v in _ff LLL commute: 

N 

h lll = (i - |n|) zid* +s™ 6 a ■ 

i=l ^ l<i<j<N 

£n a n 

Hence the ground state energy Eq (iV, a, Jl) = infspec// LLL is obtained 
from the joint spectrum of these two operators. Of particular relevance is 
the yrast curve, which is defined as the lowest eigenvalue of Ajy in the sector 
of total angular momentum L: 

Aat(L) = infspecA/v \c N =L ■ 
It is explicitly known for L < N [U[32l[20] (see also [21] for a simple proof) 

a (t\ _ ^-3/2 / lN(N - 1) for L G {0, 1} 

A 7V (L)-(2vr) / | l N{N _ 1 _ hL) ioi2 < L < N . 

Moreover, A/v(L) = for L > N(N — 1). The eigenfunctions for L = 
N(N — 1) corresponding to the eigenvalue of A^v is the bosonic Laughlin 
wave function 

\2 



n 



[Zi - 
l<i<j<N 

Little is known about Aj^(L) for N < L < N(N—1), except for numerical 
simulations for small particle number. The only rigorous result concerns the 
limit N S> 1 and L <C -/V 2 where one can show that the Gross-Pitaevskii 
approximation is exact [26]. I.e., in this regime the convex hull of An(L) is 
given by 

inf{i(/®/|*i2|/®/) : f€B, ||/|| 2 = N, (/ \zd z \ f) = L) . 
The qualitative behavior of Ajv(L) is sketched in Figure [2j 
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Figure 2. (taken from [21]) A sketch of the joint spectrum 
of Cjy and A at. The dotted line is the yrast curve, its con- 
vex hull is in red. The bold dots correspond to the possible 
ground states of ff LLL as one varies (1 — \Q)/a. The yel- 
low area on the left shows the validity regime of the Gross- 
Pitaevskii equation. For L > N(N — 1), the interaction en- 
ergy is zero. 
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